Metric perturbations in two-field inflation 
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We study the metric perturbations produced during inflation in models with two scalar fields 
evolving simultaneously. In particular, we emphasize how the large-scale curvature perturbation £ 
on fixed energy density hypersurfaces may not be conserved in general for multiple field inflation 
due to the presence of entropy as well as adiabatic fluctuations. We show that the usual method of 
solving the linearized perturbation equations is equivalent to the recently proposed analysis of Sasaki 
and Stewart in terms of the perturbed expansion along neighboring trajectories in field-space. In the 
case of a separable potential it is possible to compute in the slow-roll approximation the spectrum of 
density perturbations and gravitational waves at the end of inflation. In general there is an inequality 
between the ratio of tensor to scalar perturbations and the tilt of the gravitational wave spectrum, 
which becomes an equality when only adiabatic perturbations are possible and £ is conserved. 
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I. INTRODUCTION 

Inflation is the only known mechanism that solves the 
horizon and homogeneity problems of hot big bang cos- 
mology Q. However, the main observational prediction 
of inflationary models is the spectrum of density and 
gravitational wave perturbations they produce. Observa- 
tions of temperature anisotropies in the microwave back- 
ground, strictly speaking only provide an upper limit on 
the amplitude of such perturbations, and could in princi- 
ple be produced by some other source of inhomogeneities. 
Nonetheless, the apparently Gaussian and nearly scale- 
invariant nature of the observed perturbations are nat- 
ural properties of those produced by quantum fluctua- 
tions of the inflaton field during inflation. If inflation is 
indeed responsible for the observed anisotropies of the 
microwave background and the initial curvature pertur- 
bations from which galaxies formed, then the amplitude 
over a limited range is already constrained by observa- 
tions Q-H]. In future both the range and precision of 
these constraints promise to improve considerably and 
so it will be increasingly important to fully understand 
the predictions made by the inflationary paradigm and 
the robustness of these predictions. 

Until comparatively recently it was often stated that 
inflation predicts a scale-invariant Harrison-Zel'dovich 
spectrum of density perturbations with a negligible am- 
plitude of gravitational waves. In fact, both the tilt of 
the spectrum and the relative contribution of gravita- 
tional waves to the microwave background anisotropies 
are model-dependent quantities g]. In the conventional 
model of inflation driven by the potential energy density 
of a single slowly-rolling scalar field, the tilt and the ratio 
of tensor (gravitational wave) to scalar (density) pertur- 
bations in the microwave background can be determined 
by the slow-roll parameters which describe the slope and 
curvature of the potential as the perturbations cross out- 



side the horizon during inflation ||. Scalar curvature 
perturbations £ can be understood as originating from 
quantum fluctuations of the inflaton field that perturb 
the time it takes to end inflation ||, £ = H St — H 5<p/4>. 
Our ability to determine the linear perturbation at late 
times solely in terms of the parameters at horizon cross- 
ing depends on the constancy of £ on scales far out- 
side the horizon. The value of £ when a given comov- 
ing scale leaves the horizon during inflation can then be 
equated with that at re-entry during the radiation- or 
dust-dominated eras. 

However, most models of particle physics predict not 
only one but many coupled scalar fields and in the pres- 
ence of additional interacting fields we must re-evaluate 
these results. The conservation of Q relies on the per- 
turbations being adiabatic during inflation. In the case 
of more than one field evolving, there is the possibility 
of entropy as well as adiabatic fluctuations during infla- 
tion. We assume that all our scalar fields will eventually 
decay during reheating and only adiabatic perturbations 
remain after inflation. The validity of this assumption is 
of course very dependent upon the reheating mechanism, 
but we will leave its investigation for future work. 

In this paper we will illustrate the wider range of be- 
havior possible in multiple field inflation from the simul- 
taneous evolution of two fields. Some previous studies of 
inflation involving two fields, such as hybrid or ex- 
tended H inflation, only considered the case where one 
field evolves during inflation and the role of the second 
field is just to end inflation by a sudden phase transi- 
tion, so the single field results apply. Double inflation 
models [fiof invoke consecutive periods of inflation driven 
by two non-interacting fields. The density perturbation 
spectra produced by including an interaction term be- 
tween these fields was investigated in Ref. |ll]] . Perturba- 
tions in models that involve two interacting scalar fields 
have usually been considered in the context of Brans- 
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Dicke gravity |1 2| — 11 4|] or more general scalar-tensor theo- 
ries flBIJill , where the dilaton is expected to vary together 
with the inflaton field during inflation. However, it is the 
evolution of the second field rather than its coupling to 
the metric tensor that we wish to consider here. Only 
very recently have analytic results for general multiple- 
field inflation been presented by Sasaki and Stewart jL7) . 

We will show how to evaluate in the slow-roll approxi- 
mation the curvature perturbation at the end of inflation, 
using two alternative approaches corresponding to differ- 
ent gauge choices. In order to calculate the perturbation, 
we need to explicitly integrate along classical trajectories, 
and it turns out that it is only possible if the potential in 
separable in the different fields. We find that the pres- 
ence of entropy fluctuations modifies the usual results for 
the scalar spectra produced by inflation. In particular, 
we show that the non-conservation of £ far outside the 
horizon leads to the violation of the usual consistency re- 
lation between the ratio of tensor to scalar perturbations 
and the tilt of the gravitational wave spectrum. 



II. METRIC PERTURBATIONS 

We will consider linear perturbations about a spa- 
tially flat Friedmann-Robertson- Walker (FRW) metric 
with scale factor a(t). The most general scalar and tensor 
metric perturbations can be written as |l8|-p0t| 



ds 2 = - (I + 2A)dt 2 + 2B ti dx i dt - 
a 2 {t)\(l - 2H)Sij + 2E M - 



hij dx l dx 3 



(2.1) 



where A, B, E and TZ are scalar perturbations and hij is 
a transverse traceless tensor perturbation corresponding 
to gravitational waves. The perturbations can be decom- 
posed into Fourier modes (with comoving wavenumber k) 
which can be treated separately in the linear approxima- 
tion where they decouple. 

However not all the scalar perturbations are physical 
degrees of freedom and to remove gauge artifacts we can 
define gauge invariant quantities [[18-20 



$ = A + [B - a 2 E 



^ = TZ-H[B-a 2 E 



(2.2) 



(2.3) 



where H = a/a is the Hubble rate of expansion. Note 
that these arc equivalent to the metric perturbations A 
and TZ in the longitudinal gauge, where E = B = 0. 
Moreover, for any perturbations whose spatial part of the 
stress energy tensor is diagonal, the equations of motion 
require W = $ |20[ , so considering only linear perturba- 
tions the metric has just one scalar degree of freedom. 

The scalar TZ is the intrinsic curvature perturbation on 
hypersurfaces of fixed time t, which transforms under a 
gauge transformation, t — > t + £° as TZ — » TZ + H £° [E0| . 



It is convenient to evaluate the curvature perturbation 
on a hypersurface £ of constant energy density p,^ cor- 
responding to the choice of gauge £° = Sp/p pi] . This 
gives TZy: = TZ + H Sp/p. The intrinsic three-curvature 
on this surface is = 4V 2 Kj. Note that if we start 
from a flat hypersurface (TZ = 0), we can interpret TZs as 
being due to a perturbation in the logarithm of the scale 
factor (or number of e-foldings N) on that hypersurface 



TZv = SN = HSt , 



(2.4) 



where St = Sp/p. In particular for inflation with a single 
field we can write St = S<p/<p, which gives the familiar 
result for the origin of density perturbations Q. 
We will define a quantity 



v 2 $ 

3H 



H 2 

$ s- («]> - II [ «J> 



(2.5) 



written in terms of the gauge invariant metric potential, 
which coincides with TZs on scales far outside the horizon 
(k <C aH). The time-dependence of ( on these super- 
horizon scales is given by 



(2.6) 



where Sp and Sp are the perturbations in the pressure 
and energy densit y o n spatially flat hypersurfaces. As 
we can see, Eq. ( |2.6| ) vanishes for adiabatic perturba- 
tions. This is the case for inflation with a single field, so 
we can evaluate the curvature perturbation at late times 
by equating it with that at horizon crossing. In fact, £ 
continues to be conserved on super-horizon scales during 
the radiation- and dust-dominated eras, and therefore, 
we can directly equate the curvature perturbation when 
it left the horizon during i nfla tion with that at re-entry. 

The expression in Eq. (2.5) allows one to evaluate £ 
once we know To that end, one can integrate the cou- 
pled perturbation equations for $ and the scalar fields 
in a particular gauge, see for instance Ref. This 
allows us to give an expression for £ at the end of infla- 
tion in terms of the field fluctuations at horizon crossing. 
Alternatively, one could find £ by evaluating the pertur- 
bation in the number of e-foldings, SN, given in Eq. 



C ~ TZz = SN , 



(2.7) 



as proposed recently by Sasaki and Stewart |Tj. On 
spatially flat hypersurfaces the perturbed fields on super- 
horizon scales effectively obey the same equations of mo- 
tion as the homogeneous background fields, see Eq.(|3.2|). 
Thus SN is the difference in the number of e-foldings 
along neighboring classical trajectories in phase-space. In 



*In single-field inflation, this corresponds to a comoving 
hypersurface. 
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order to evaluate 5N at the end of inflation in terms of 
field fluctuations at horizon crossing, we have to integrate 
the background equations of motion, not only along a sin- 
gle trajectory but also along the perturbed trajectories. 
We will show that this does indeed yield the same results 
as those obtained by directly solving the evolution of the 
metric perturbations in the longitudinal gauge. The dif- 
ference between the two approaches is just a choice of 
gauge. 

In practice, whichever method one adopts, one needs 
to know not only the initial perturbation but also its in- 
tegrated effect along the subsequent trajectory. In the 
following sections we shall show how the perturbations 
can be evaluated in specific models. However it is im- 
portant to remember that £ only remains a conserved 
quantity thereafter if we can treat the end of inflation 
as a transition at fixed energy density and the subse- 
quent evolution is adiabatic. If we wish to match metric 
perturbations at the end of inflation, across a hypersur- 
face of fixed energy density, we must match £ on large 
scales |2l| . While we expect this to be true in single field 
models of inflation, it is a much more complicated issue 
in two-field inflation, since it depends upon the dynam- 
ics of reheating. Such an involved issue deserves further 
attention and is beyond the scope of this paper. 

The attractive feature of Sasaki and Stewart's ap- 



proach is that Eq. (2.7) is a purely geometrical result, 
independent of the matter content (subject only to the 
condition W = $) and could in principle also be applied 
to calculate £ on surfaces long after the end of inflation. 
However evaluating 8N in the radiation-dominated era 
once again requires a quantitative understanding of re- 
heating along different trajectories. In what follows we 
will restrict ourselves to a calculation of the curvature 
perturbation at the end of inflation. 



III. TWO-FIELD MODELS OF INFLATION 

In this section we will consider a model with two scalar 
fields, described by the action 



S = 



for instance, in the Einstein frame[j] of general scalar- 
tensor gravity theories Jl6|], in which case W((f>,a) — 



-4h(4>) 



V(a) 



The field equations for the fields a and in a spatially 
flat FRW metric are then 



o + 'iHa 



Z 2H W dW +2h , w ^ & 
do 



dW 

3H(f>=-^--ti(4>)e 



-2h(4>) ,5-2 



(3.2) 



and the Hamiltonian constraint is 

..2 



H 

~6~ 



4? + e ~ 2hW & 2 +2W(<j),o) 



(3.3) 



The condition for inflation to occur \H\ < H 2 is thus, 
see Eqs. Q and (jO), 



" 2 + e- 2h & 2 < W(<f>,a) . 



(3.4) 



In calculating perturbations on a comoving scale k, we 
shall see that an important quantity is the number of e- 
folds from the end of inflation when that scale crossed 
outside the horizon (fc = a*i?*), 



N = - 



H dt. 



(3.5) 



Our present horizon crossed outside the Hubble scale 
about 50 to 60 e-foldings before the end of inflation. The 
precise number depends logarithmically on the energy 
scale during inflation and the efficiency of reheating, and 
so is weakly model-dependent. 

We can derive an exact expression for the time depen- 
dence of £ using the linearly perturbed field equations 
@ 



(3.6) 



where 



-e-2 h Wg ab a, a a, b -W(<t>,<T) 



(3.1) 



where R is the usual Ricci curvature scalar and k 2 = 
8irG. If h = then the fields have standard kinetic terms, 
but we have also allowed for the possibility that the a ki- 
netic term has a ^-dependent pre-factor as would come 
from a conformal transformation of a theory with a non- 
minimally coupled <fi field. Such an action might arise, 



2 dt y 02 + e - 2h & 2 J \(t> 2 + e- 2h o 2 J 

(3.7) 

Note that although 5<j) and 5a are gau ge dependent quan- 
tities, their combination in Eq. (3.6) is gauge invariant. 



^The original Jordan frame, in which a is minimally coupled, 
is related to the Einstein frame used here via the conformal 
transformation g a b = e 2h ^ g a b- 
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If only one field is evolving (<r = or <j> — 0), we see that 
T = and ( is conserved on large scales (k <C aH), 
thus recovering the well known result p2|. However, 
if T 7^ and both fields are evolving, £ only remains 
constant for perturbations along the classical trajectory 
(Sa/Scj) = &/<p). The first term on the right-hand-side 



of Eq. (3.7) is present whenever two fields are evolving. 
The second term is due solely to the presence of the non- 
standard kinetic term for a and vanishes when h = con- 
stant. It represents the frictional damping of the a field 
by h{4>). 



A. Slow-roll trajectories 

To make further progress we will work in the slow- 
roll approximation in both scalar fields. In principle this 
is not necessary for inflation to occur: one of the fields 
might roll quickly to the minimum of its potential and 
then the problem reduces to single field inflation. Models 
of hybrid inflation or other models of first-order infla- 
tion ^|J|] provide examples where more than one field 
is present but only one field slow-rolls during inflation. 
However, here we wish to consider the case in which both 
fields slow-roll. The slow-roll ap prox imation amounts to 
reducing the full field equations (3.2) to first-order equa- 
tions, 



3H 2 



K 2 W{(t>,(j) 



2h(d>) 



d\nW H 



da k 2 
d\nW H 



d<j) k 2 



(3.8) 
(3.9) 

(3.10) 



This approximation reduces the effective four-dimen- 
sional phase-space during inflation down to the two- 
dimensional field space {0, a}. 

This leads to a crucial difference between single-field 
inflation and inflation with two or more fields. In 
single-field inflation the slow-roll solution forms a one- 
dimensional phase-space, i.e., there is a unique trajectory. 
The end of inflation corresponds to a fixed value of the 
field and any comoving scale which crossed the horizon N 
e-foldings before the end of inflation also corresponds to 
a unique value which may be calculated from Eq. (3.5). 



In two-field inflation the slow-roll approximation leaves a 
two-dimensional phase-space. Classical trajectories dur- 
ing inflation in this field space correspond to lines which 
are no longer unique. In particular, the end of inflation 
will in general be described by a one-dimensional line in 
this phase-space rather than a single point, as will the 
locus of a given number of e-foldings from the end of 
inflation. 

To calculate £ using the formalism of Sasaki and Stew- 
art we need to know how the dependence of the number 
of e-foldings from the end of inflation upon the pertur- 
bations in the fields (j> and a. In the single field case 



this can only amount to a perturbation along the clas- 
sical trajectory, St — 5<f>/(f>, due to the equations of mo- 
tion 0. However, in the two-dimensional field space, 
perturbations will in general move the fields onto a dif- 
ferent trajectory with a different end point of inflation, 
except when the perturbations happen to be adiabatic. 

Therefore, it may no longer be sufficient to use the 
familiar relations along a given trajectory 7, 



ON 



H 



ON 
~da~ 



H 

a 



(3.11) 



in order to work out the change in the number of e- 
foldings due to perturbations in t he fi elds. To evaluate £ 
at the end of inflation using Eq. ( |2.7D , we have to allow 
for variations away from the classical trajectory. 
Fortunately, in the case of a separable potential, 



W{M = U{4>)V{a), 



(3.12) 



we can label the slow-roll trajectories by an integral of 
motion, 



C 



V da 



Ud<t> 



(3.13) 



which allows us to parametrize motion off the classical 
trajectory. We can then substitute the slow-roll equation 
of motion for <f> into Eq. (3.E) to obtain the number of e- 
foldings for a given value of </)* along a trajectory labeled 
by C, 



U, 



U'(cf>) 



(3.14) 



This depends on the values of both (/>* and er* through 
the dependence of <f> e upon C ((/>*, it*). 

One should, of course, be cautious about evaluating 
the perturbation by integrating along a slow-roll trajec- 
tory that is only an approximation to the full equations 
of motion. In the single-field case the conservation of £ 
is an exact result and does not rely on the slow-roll ap- 
proximation. In the examples presented in this paper we 
have checked by numerical solutions that the slow-roll 
results remain a good approximation for a wide range of 
parameters right up until the end of inflation. 



B. Scalar perturbations 

Perturbations during inflation are expected to arise 
from the vacuum fluctuations of the fields which are 
stretched by the inflationary expansion up to super- 
horizon scales. In the following we refer to the gauge in- 
variant scalar field perturbations p0[ , or equivalently the 



'One can verify from the slow-roll equations of motion that 
C ~ along any classical trajectory 7. 
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scalar field perturbations in the longitudinal gauge. For 
large values of k aH we can neglect the potential terms 
in the perturbed field equations for the scalar fields and 
they reduce to those of massless fields. Thus, to lowest 
order in the slow-roll parameters, the expectation values 
of the perturbations as they cross outside the Hubble ra- 
dius (k ~ a»ii*) are given by Gaussian random variables 
with e~ 2h * (|,5a* | 2 ) = Hi /2k 3 , (|^| 2 > = H 2 /2k 3 , 
where k is the comoving wavenumber. Note th at, while 
the field 4> has a standard kinetic term in Eq. (3.1), the 
a field does not and so the ^-dependent prefactor must 
be included in the expectation value acquired at horizon 
crossing We shall denote the spectrum of a quan- 

tity A by V A (k) = (47rfc 3 /(27r) 3 ) (\A\ 2 ) , as defined in @. 
Thus we have 



V &4> 

Vsa 



2^ 

„2/i„ 



/ / 

2^ 



(3.15) 
(3.16) 



For slowly varying, long-wavelength (fc <C aH) modes, 
to lowest order in the slow-roll parameters, the pertur- 
bation equations can be integrated to give |1(| , 



8<j>: 



1 U'{4>) 



Sa^ -1 



U(<P) 
V(a) 



(Q 2 



'2h 



Qi 



(3.17) 



(3.18) 



where Q\ and Qi are constants of integration. It will be 
convenient to define a new constant Q3 = Q2 + e 2h *Qi, 
so that Qi and Q3 are independent Gaussian random 
variables whose values, for a given Fourier mode, are de- 
termined by the amplitude of 5a* and Scj)* at horizon- 
crossing. Thus they have expectation values 



=-( 


'HA 


H \ 




k 2 e 2h * 


(H* 


24 


\2ir 



(3.19) 



(3.20) 



During slow-roll, in the long wavelength limit, the cur- 
vature perturbation on hypersurfaces of constant energy 
density can be written as GaM 



-2h- 



aSa 



[e + ( 



e 2h -e 2h ^ 



e<r] Ql + t a Qz 



(3.21) 



+ e 2h e a 

where we have extended the usual definition of the slow- 
roll parameter e for a single field p5| to 

2 -, /TTll IS \ 2 

, £0 



1 (V'(a) 



2k 2 V V{a) 



1 (u'{<t>)X 

2k 2 \ U{4>) ) 



(3.22) 



If either of the scalar fields is fixed, t a or identically 
zero, then we recover the single field results where £ is 
constant and equal to Qi or Q Sl respectively. We see 
from Eq. (3.6) that C continues to evolve after horizon 
crossing if both fields are evolving and T / 0. In the 
slow-roll approximation we have 



c 



£{h,-h) 



-2h, 



Q3-Q1) T. 



(3.23) 



However, if h grows significantly during inflation, then 
the change in £ may be small, as can happen in the case 
of scalar-tensor gravity theories j^O]. On the other 
hand, for minimally coupled scalar fields, where h = 0, 
we have A£ ~ (Q3 — Qi) J Tdt which is typically of order 

C*- 

The spectrum of density perturb ation s at the end of 
inflation V^k) computed from Eq. ( 3.21 ) is, 



C 2 I 2vr 



+ (e 



2/1, 



,2h» 



c2h. 



a 2h. 



1 



(3.24) 



where we have used Eqs. (|3.19| ) and (|3.20[) . This is the 
result found in Ref. Gq|. 



We will now show that Eq. ( 3.24 ) can also be derived 
in the framework of Sasaki and Stewart |l7j in terms 
of the change in the number of e-foldings, 8N , as given 
in Eq. fl2.7p. We wish to evaluate £ e , the curvature 
perturbation on a hypersurface of fixed energy density, 
U(<p e )V(a e ) = constant, near the end of inflation. The 
values of the fields cj) e and a e at the end of inflation on 
a given classical trajectory will be a fun ction of the con- 
served quantity, C, defined in Eq. ( 3.13| ). Therefore, 



dC 



da e 
dC 



0. 



(3.25) 



Differentiating Eq. ( 3.13 ) with respect to the trajectory 
C, and using Eq. (3.25), we find 



da e 
dC 



dC 



„2h e 



„2h e 



(3.26) 



We can now ev aluat e the dependence of the number of 
e-folds TV, Eq. (3.14), on the initial values <j)„ and er*, 



dN 



d(f>* 



dC 



dC 



+ ^ — da * 
oa* 



(3.27) 
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Using Eqs. (|3.26[) and (|3.13|) we find 



dN 


2 


n 2 




- (e 2,l = - e 2 ' 1 * 


K 






Of* 










dN 


2 


K 2 


C 


2 




da* 




2e* 







(3.28) 



These are the expressions required to evaluate Vf us- 



ing Eq. ( 2.7) with 
and (3.16), 



and 8a* given by Eqs. (3.15) 



m 




dN 


2 

+ £ 2K 


ON 


2" 






dcj>* 




da* 





(3.29) 



This expr ession exactly coincides with that obtained in 
Eq. (3.24). Note that this result applies to both scalar- 
tensor theories (with U (4>) = cxp{— 4h((j))}, see Ref. @), 
and for minimally coupled (h((/)) — 0) two-field inflation 
with a separable potential in general relativity. 



C. Tensor perturbations 

In addition to the scalar curvature perturbations that 
give rise to density perturbations, ten sor or gravitational 
wave perturbations [hij in Eq. (2.1)] can also be gen- 
erated from quantum fluctuations during inflation p6[ . 
These do not couple to the matter content and are deter- 
mined only by the dynamics of the background metric, 
so the standard results for the evolution of tensor per- 
turbations of the metric remain valid. The two indepen- 
dent polarizations evolve like minimally coupled massless 
fields with a spectrum [EOLbI 



(3.30) 



Note that the tilt of the gravitational wave spectrum, 
n g = dlnVg/dhik, is given by 



(3.31) 



Noting the condition for inflation given in Eq. (3.4) we 
see that the definitive test for inflation is the presence of 
a gravitational wave spectrum remains —2<n g < and 
is unaltered by the presence of more than one scalar field. 
The actual measurement of this slope will be exceed- 
ingly difficult. Tensor perturbations do not contribute 
to structure formation and in many inflationary models 
the observable effect of gravitational waves is completely 
negligible 

Gravitational wave perturbations can contribute to the 
microwave background anisotropies only on the largest 
scales (scales larger than the Hubble scale at last- 
scattering, corresponding to about 1° on the sky). Their 



contribution relative to scalar curvature perturbations is 
given by the ratio [|| 



R -lv~ c 



(3.32) 



The rapid decay of the gravitational wave anisotropies 
on smaller scales is their most distinctive signature. If 
we define 



G=-^+e 



2h 



(3.33) 



then we can write the ratio of tensor to scalar contribu- 
tions as 



R~12el 



G 2 (G* 



,2h* 



" (Ge 



a 2h, 



{G* - e 2h ') 



(3.34) 



Together with Eq. ( 3.31 ) one can show, see Fig. 1, that 

R < G\n g \ . (3.35) 

This result was found by Polarski and Starobinsky in 
their model of double inflation |27]] and given in Ref. JlJ] 
for general multiple-field inflation. It generalizes the 
usual consistency relation between R and n g in single 
field inflation, where R ~ 6|n s |, see Ref. pifl . It shows 
that in principle one could tell from the spectra of met- 
ric perturbations whether more than one field is evolving 
during inflation. 



The inequality of Eq. (3.35) beco mes a n equality only 
when G* = G e . Note that C in Eq. (3.21) can be written 



as 



c 



(G-e 2h *)Q 1 + Q 3 
G 



(3.36) 



and therefore the inequality is saturated whenever £ re- 
main s co nstant after horizon crossing. As shown in 
Eq. ( [2.6| ) this only occurs when the perturbations are 
adiabatic. This is, of course, true when only one field is 
evolving but we can also find a class of two-field models 
in which G is constant in the slow-roll approximation. 

Even perturbations off the classical trajectory are 
forced to be adiabatic (dp/ Sp = p/p) when the pressure p 
is a function solely of the density. In our two-field model, 
the pressure and energy density are given by, 



1 

P= 2 



P= 2 



-2h ■ 2\ 

e a 



W{4>,a), 



+ e- 2h a l ) + W{<p, a), 



(3.37) 



(3.38) 



and so p = p(p) implies that the kinetic energy density 
is a function of the potential. We will now show that in 
the slow-roll approximation for a separable potential this 
leads to G being conserved. Along the classical trajecto- 
ries, 



G 



2h'e 



I „2h 



(3.39) 
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If the kinetic energy is to be a function solely of the 
potential energy density, we require + e 2h e a — f(W), 
which implies that 



G = 



dlnW 



& dlnW 



df 



da 



d\nW 



(3.40) 



which from the slow-roll equations of motion must be 
zero. Thus we see that both the relation R = 6\n g \ 
and the constancy of £ are a consequence of the scalar 
perturbations being adiabatic rather than the number of 
fields present. 

For example, consider the class of generalized Brans- 
Dicke theories p9j , |l4| ] where 2h{<f>) = jutfi and U{4>) — 
e -0K<f> m j.jjg jr; ms t em frame with a polynomial inflaton 
potential, V(a) — Xa 2n /2n. In this case, we have 



£0 + e 2h e B 



2 



2n 2 



(3.41) 



We see that for n — we have f(W) = /3 2 /2 + 

2r? I k 2 {2nW I X) 1 1 n . As discussed above, this ensures 
that all perturbations are adiabatic and thus C is con- 
served. We could also see this by evaluating G, which 
in this case becomes G = (3jC, where C is the con- 
serve d qu antity along the classical trajector y, gi ven in 
Eq. ( 3.13| ). Thus C is constant, see Eq. (3.36), and 
R = 6 \n g \. 

In particular, we find that £ remains a constant after 
horizon crossing in Brans-Dicke gravity, where (3 = 27, 
for a quartic potential V = Aer 4 /4. We will study this 
case in more detail in the next subsection. 



D. Brans-Dicke case 

We will study here a very simple case which we can 
solve completely. This is Brans-Dicke case with a quar- 
tic potential for the inflaton field, h{4>) — aK<j), £/(</>) = 
e -4afc<^ anc j V(a) = Act 4 /4. The constant a characterizes 
the relative coupling of scalar and tensor fields to mat- 
ter, and is related to the usual Brans-Dicke parameter by 
2a 2 = 1/ (2uj + 3) . As shown in the preceding section 
this belongs to the sub-class of two-field models for which 
C is in fact conserved on scales outside the horizon in the 
slow-roll approxi matio n. 

Following Eq. ( 3.13 ), each classical trajectory can be 
parametrized by 



G = 8a C = — + e 

Err 



2/i 



(3.42) 



where = 8a 2 and e 



rameters defined in Eq. ( 3.22 ). 

The spectrum of curvature perturbations is then 



k a are the slow- roll pa- 



V',- 



k 2 fH*\ 2 1 



2lT 



Gt 



k 2 /i3* x 2 



2tt 



(3.43) 



This is an extremely simple and compact formula which 
makes it possible to compute the spectral tilt of the scalar 
perturbations, n s — 1 = dlnV^/dlnk, 

n s -l = ~2 H - 3e 2h * e* . (3.44) 
The ratio of tensor to scalar perturbations is given by 



R~ 12 Get ~ 12(e d) + e 2h 'e 



(3.45) 



Depending on e 2h e a at horizon crossing, we might have 
significant gravitational wave contributions. As shown 
earlier, the tensor perturbations' spectral index is given 
by n g ~ —R/Q, the same relation as in the single field 
case. 



E. Minimally coupled case 

If we restrict our attention to fields which are mini- 
mally coupled, i.e., h = 0, with a separable potential, 
then there is no frictional damping of the a field and 
the results simplify compared with the general case. The 
spectrum of curv ature perturbations at the end of infla- 
tion in Eq. (3.24) becomes 



n 2 fH* 
2tt 



1 



1 



(3.46) 

and the ratio between the tensor and scalar contributions 
to the microwave background anisotropies on large scales 
is given by 



R~ 12 



{e% + e%) 2 e% e* 



(el) 2 el + [e%) 2 e* 



(3.47) 



Only if either or e„ are zero, or if both are constants 
(corresponding to exponential potentials) , are all pertur- 
bations adiabatic, Q conserved and R = 6 \n g \. 

More generally, the perturbations may still become ef- 
fectively adiabatic by the end of inflation if the evolu- 
tion has become essentially one-dimensional along, say, 
the a direction in field-space. This requires not only 
e ff ^ e d>' a lso C on a given scale to be due to per- 



turbations in the a field, which from Eq. (3.21) requires 

( % > e % (Q1/Q3) - (-% \l e *cr/ e %- In this casc C has be- 
come constant by the end of inflation and, just as in the 
single field case, we expect it to remain conserved until 
that scale re-enters the horizon during radiation- or dust- 
domination. Nonetheless we may still see evidence of the 
second field evolving at horizon crossing during inflation 
due to the decrease in the Hubble rate, H*, caused by 
the slow-roll of tj> in addition to that due to the evolution 
of a. The tilt of the scalar perturbation spectrum in this 
limit is 
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l~-6e; + 27j*-2ea 



where we have introduced a further slow-roll parameter 
r\ a = V" /k 2 V which describes the curvature of the po- 
tential along a. Only as — * do we recover the familiar 
single- field result |^8|,|| . The tilt of the spectrum of ten- 
sor perturbations is n g = — 2(e* + et). Thus the slope of 
the second field tends to decrease both n s and n g . How- 
ever the ratio between tensor and scalar modes is given 
by R ~ 12e* and so is a function solely of the slope of the 
potential along a at horizon crossing, as in the single field 
case. Thus the evolution of the <f> leads to the violation 
of the single-field consistency relation, R ~ —6n g . 



IV. CONCLUSIONS 

The presence of more than one field evolving during 
inflation has important consequences for some of the fa- 
miliar results quoted in inflationary cosmology. In the 
slow-roll approximation, n fields lead to an n-dimensional 
phase-space. Thus, there is no longer a single classical 
trajectory that leads to the end of inflation and different 
initial conditions may lead to different end-points. Quan- 
tum fluctuations lead to perturbations not only along 
the trajectory but also onto neighboring trajectories. We 
have entropy as well as adiabatic perturbations and, in 
this case, the curvature perturbation £ is no longer con- 
served on super-horizon scales. Therefore the perturba- 
tion on a given scale is no longer determined solely by 
quantities at horizon crossing, but also depends upon the 
subsequent evolution. In order to evaluate £ at the end 
of inflation we must be able to integrate the perturba- 
tion along the trajectory. We have shown how this may 
be performed for a separable potential in the slow-roll 
approximation. 

In single-field inflation the ratio between the contri- 
bution of tensor to scalar perturbations of the microwave 
background on large scales can be related to the tilt of the 
gravitational wave spectrum, R = —6n g |2^||. We have 
shown that this is a consequence of the scalar perturba- 
tions being adiabatic and can also occur, within the slow- 
roll approximation, in a limited class of two-field models. 
More generally, in two-field models with a separable po- 
tential the ratio obeys the inequality R < Q\n g \, wh ere n a 



(3.48) evolution after inflation and in particular to the dynam- 
ics of reheating. If the end of inflation corresponds to a 
fixed energy density then £ will be conserved across this 
hypersurface. This is the case when there is a unique 
end-point and only one field is evolving at the end of in- 
flation. However in the presence of two or more fields 
this need no longer be true and requires a more detailed 
study. 

An important question remains as to whether, even 
if a given theory has many scalar fields, one should ex- 
pect to see evidence in the perturbation spectra of more 
than one field evolving. The evolution during the final 
60 e-foldings of inflation depends upon the initial condi- 
tions, and this may lead to only one field evolving at late 
times. For instance, in the case of chaotic inflation in 
some scalar-tensor gravity theories, the stochastic evo- 
lution may fix the value of the gravitational coupling, 
effectively leading to single field inflation 30 1 . To answer 
the question of initial conditions, one needs to under- 
stand the stochastic evolution of the coupled fields at 
early times |l|M||. 



is given in terms of the slow-roll parameters in Eq. ( 3.31 ) . 
While in principle this is a distinctive prediction of infla- 
tion with more than one scalar field, it is very difficult 
to observe. By contrast the tilt of the scalar spectrum is 
already constrained by observations but does not give a 
model-independent test of multiple-field inflation. 

While we have been careful to calculate the evolution of 
the curvature perturbation during inflation, we have not 
attempted to go beyond the end of inflation. Since the 
observable quantity is the amplitude of scalar perturba- 
tions at re-entry during the radiation- or dust-dominated 
eras, the interpretation of our results is sensitive to the 
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FIGURE CAPTIONS 



Fig. L Plot of the ratio R/6\n 



and (3.34) as a function of the parameters a 



1 and b = G* e 
imum of one for a 
the minimally coupled case (h 
so both a and b must be non-negative 



from Eqs. (3.31) 
G e e^- 

1 = e^/e*. This ratio reaches a max- 
b, corresponding to G e = G*. In 
0), we have a = e%/e%, 



<) 



